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Experimentally, ultracold (T ^ 1 K) neutral plasmas could be realized only recently 
by photoionizing a cloud of ultracold atoms [T1I21IS]- From a plasma physics perspective, 
they are very appealing since they might provide laboratory realizations of so-called 
strongly coupled plasmas, where the Coulomb interaction (inversely proportional to the 
mean inter-particle spacing a) dominates the random thermal motion (proportional to 
the temperature T). As a consequence, interesting ordering phenomena such as Coulomb 
crystallization can be observed. In order to be in this strong-coupling regime, a system 
has to be either very dense (large Coulomb interaction) or very cold (little thermal 
energy). The first scenario is realized in some astrophysical contexts, which are difficult 
to access experimentally. Hence, the possibility to realize a strongly coupled plasma 
in the laboratory by going to ultralow temperatures offered exciting prospects within 
the plasma physics community. In [T] , an ultracold neutral plasma has been created by 
photoionizing a gas of laser-cooled Xe atoms. The extremely low temperatures of the 
Xe gas suggested that the system was well within the strong-coupling regime, the so- 
called Coulomb coupling parameter T = e^/ {ak-eT), i.e. the ratio of potential to thermal 
energy, being of the order of 10 for the electrons and even up to 30000 for the ions|. 

However, it was realized quickly [H E] that despite the low temperature of the 
Xe gas a strongly coupled plasma could not be created in this way. Due to the fact that 
the neutral Xe atoms initially present interact only very weakly, the plasma is created 
in a completely uncorrelated non-equilibrium state. Hence, the subsequent conversion 
of potential into kinetic energy rapidly heats both the electron and ion subsystem, 
suppressing the development of substantial correlations [H Ej. This effect has been 
termed "disorder-induced heating" [H] or "correlation-induced heating" [S], expressing 
the fact that it is the development of spatial correlations during the relaxation of the 
plasma towards an equilibrium state which leads to an increase in temperature. Put 
another way, in the initial non-equilibrium state the potential energy is higher than in the 
corresponding equilibrium state of a plasma of the same average density. Consequently, 
relaxation towards thermodynamic equilibrium will reduce the potential energy in the 
system and increase the kinetic energy, i.e. the temperature§. Additionally, the electrons 
are heated by three-body recombination as well as continuum lowering. Suppression of 
electron heating is thus much more difficult than simply reducing the correlation heating. 
At present, there is a proposal of mixing the plasma with a gas of Rydberg atoms, so 
that coUisional ionization of these Rydberg atoms removes energy from the electron 
component of the plasma jH]. 

For the creation of a strongly coupled ionic component, different procedures have 

X The large mass difference between electrons and ions implies a very long timescale for energy transfer 
between the two subsystems. Hence, for practical purposes it is justified to treat them as (almost) 
independent and assign different temperatures to both subsystems. 

§ In view of the fact that the initial state of the plasma is far from equilibrium, it is not really justified to 
associate a temperature with it, the properties of the system are not those of a state in thermodynamical 
equilibrium with the same kinetic energy. However, before the photoionization occurs the atomic gas 
has a well-defined temperature, and it is the temperature of these atoms that can be compared to that 
of the ions after equilibration. 
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been suggested to avoid these heating effects, namely (i) using fermionic atoms cooled 
below the Fermi temperature in the initial state, so that the Fermi hole around each 
atom prevents the occurence of small interatomic distances [Tj; (n) an intermediate 
step of exciting atoms into high Rydberg states, so that the interatomic spacing is at 
least twice the radius of the corresponding Rydberg state [Hj; and {iii) the continuous 
laser-cooling of the plasma ions after their initial creation, so that the correlation 
heating is counterbalanced by the external cooling ^21 So far, none of these 

procedures has been realized experimentally, however it has been shown theoretically 
that continuous laser-cooling can indeed lead to a strongly coupled plasma, connected 
with crystallization of the expanding system In this case, the correlation heating 
effect is still present and is offset by an additional external cooling mechanism. Proposals 
(z) and (ii), on the other hand, aim at suppressing the heating directly by avoiding 
the small interparticle distances leading to a large potential energy, i.e. by introducing 
(spatial) correlation in the initial state. As has already been noted in |B] , this could also 
be achieved using optical lattices to arrange the atoms. 

At this point, the question arises naturally whether the correlation heating could 
not only be avoided, but whether it can actually be reversed and turned into cooling. As 
argued above, correlation heating describes the fact that the uncorrelated ions created in 
the photoionization process have a higher potential energy than in the thermodynamical 
equilibrium state. Hence, equilibration will reduce the potential energy and increase the 
kinetic energy of the system, making the ions hotter than the atoms in the initial gas 
state. On the other hand, it is conceivable that the atomic gas might be prepared in an 
"overcorrelated" state, where the potential energy is lower than its equilibrium value. 
In this case, equilibration must increase the potential energy on the cost of the kinetic 
energy, i.e. the temperature must decrease. This effect of "correlation cooling" has been 
described before in |^ and, based on a kinetic approach, in more general terms in |13j . 
In [12], the overcorrelated initial state was achieved by suddenly decreasing the Debye 
screening length, i.e. the temperature of the electrons, rather than by inducing spatial 
order. It was shown that in this way, a reduction in temperature of the order of a few 
percent does indeed occur for the temperatures and densities considered. 

In the present letter, we explore an alternative scheme for achieving correlation 
cooling, namely the use of an optical lattice to induce spatial correlations in the initial 
state as suggested in p. We describe, both analytically and by molecular dynamics 
simulations, the situation of a fully ordered initial state where the atoms have been 
arranged in a perfect bcc-type lattice structure using an optical lattice. In reality, 
their spatial distribution will be broadened, classically due to their non-vanishing 
kinetic energy corresponding to the initial-state temperature, quantum mechanically 
corresponding to the eigenstates of the harmonic oscillator potential generated by the 
optical lattice. This broadening is neglected in our numerical simulations, however, 
its influence will be discussed in the analytical considerations below. For the sake of 
numerical simplicity, our initial state corresponds to a sphere cut out from a bcc lattice 
rather than a more realistic cylindrical shape, arguing that edge effects should be of 
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minor importance if the number of atoms in the system is large enough. The system is 
propagated using the hybrid molecular dynamics method described in ^1]. Briefly, an 
adiabatic approximation is employed for the electrons, and their distribution is obtained 
from the steady-state King distribution [T3j , widely used in the study of globular clusters. 
The ions are then propagated using the electronic mean-field obtained from the steady- 
state distribution and by fully taking into account the ion-ion interactions. In the 
experiments under consideration, quasineutral plasmas with Fg as low as 4 ■ 10~^ j2] 
have been realized, making Debye screening effects negligible, while ion numbers iVj 
up to 10^ have been obtained For our simulations, such a large particle number 
would lead to a prohibitively large numerical effort. However, the typical timescale 
for the plasma expansion is texp ~ V^N^^^Uq^, where c^o = \J /Ma^ is the Einstein 
frequency in the Wigner-Seitz model of a Coulomb solid and M is the ion mass. On the 
other hand, the relaxation of the ions takes place on a timescale of trei ~ ^o^- Hence, 
the same macroscopic behavior can be obtained by decreasing iVj and increasing Fg, 
keeping the ratio tgxp/^rei fixed. Therefore, for our simulations we chose A^i = 10^ and 
an initial electron energy corresponding to Fg = 0.5, while the bare Coulomb potential 
was used for the ion-ion interaction. 

Figure ^a) shows the time evolution of the average ion kinetic energy in units of the 
initial-state temperature for a typical realization. Since the initial lattice configuration 
produces an effective oscillator potential for each ion, the relaxation process is connected 
with transient oscillations of the kinetic energy which are damped out due to ion-ion 
collisions. Additionally, the average kinetic energy increases at later times due to a 
radial ion acceleration caused by the thermal electron pressure, which is proportional to 
the density gradient. Therefore the slow expansion of the unconfined plasma leads to a 
softening of the initially sharp plasma edge and a modification of the velocity distribution 
in the outer region of the plasma only. Since this distribution is close to a Maxwellian 
in the central region of the plasma, the average kinetic energy is directly related to 
the ionic temperature T, while it becomes more and more dominated by the expansion 
energy towards the edge of the plasma. In order to highlight these edge effects, we have 
calculated the average ion kinetic energy from the velocity distribution in the inner 
region only, taken to be a sphere with half of the initial plasma radius. The resulting 
curve is shown as the dashed line in figure^a). Clearly, the expansion hardly affects the 
kinetic energy in the inner region of the plasma, leaving only the oscillations which are 
damped out as the system approaches thermodynamic equilibrium. As discussed above, 
the relaxation of the ions towards equilibrium takes place on a timescale of oJq^, while 
the typical timescale for the plasma expansion is \/TIn1^^ujq'^ . Hence, a sufficiently 
large number of ions allows for an almost unperturbed relaxation in the inner plasma 
region so that a temperature can be defined in a meaningful way. (For a strontium 
plasma with a typical density of 5 ■ 10^cm~3 used, e.g., in [12], uJq^ ~ 0.2 /is, while 
tgxp is about 30 times larger for Fg = 4 ■ 10~^ and A^j = 10^.) Using the absorption 
imaging techniques described in [12] it should be possible experimentally to probe this 
inner region only. Therefore, we chose to define the final temperature as given by the 
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Figure 1. (a) Time evolution of tlic ionic kinetic energy, for the whole plasma cloud 
(solid) and for the central region only (dashed), (b) Ratio of final to initial temperature 
as a function of the density-scaled inverse initial ion temperature Fj, from the numerical 
simulations (full circles) and from equation © (solid line). The inset shows the final 
Coulomb coupling parameter. 



inner plasma region and to neglect the finite-size effects associated with the boundary 
region of the plasma. The resulting Tf was estimated from the center of the remaining 
oscillations and is shown in figure ^b) as a function of Fi by the full circles. It should 
be noted that the value of Fi has no meaning as a measure of nonideality of the initial 
plasma state, since the system may be far from equilibrium. Therefore, in the present 
case it should simply be viewed as a density-scaled inverse temperature rather than a 
Coulomb coupling parameter. As can be seen in the figure, the ratio of final to initial 
temperature reaches a constant value of about 0.5 as the initial temperature goes to zero. 
Hence, cooling can indeed be achieved on a moderate level of a temperature reduction 
of about 50 percent. 

For large particle numbers and if Fg is small enough such that three-body 
recombination and Debye screening can be neglected, quasineutrality leads to a 
homogeneous electron background in the central plasma region. Thus, we can use the 
one-component plasma model to gain more insight into the cooling process. This model 
has been studied extensively over the last decades [ini IIZ|, which led to a detailed 
understanding of the equation of state over a wide parameter range. In the classical 
regime, the energy per particle of the final equilibrium state is given by 

f/f = ^k^Ti + Wi = i^ + w,)k^T, , (1) 

where the excess energy W accounts for nonideal effects due to the inter-particle 
Coulomb interaction. For this quantity, several fit formulae can be found in the 
literature, which have been obtained by Monte Carlo and molecular dynamics methods 
[TBI 113 HE] • In the fluid phase we use the interpolation formula from [12] 

where Ai = -0.9052, A2 = 0.6322 and A3 = -^3/2 - Al/^/A^. Equation ^ satisfies 
the known small-F limit [20] and accurately fits the high-F values obtained by numerical 
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calculations |18|. In the solid phase (F > 174), the excess energy can be described by 
an expansion in powers of 1 /F jJTj 

(3) 

f If ■ " ■ ■ 

where a = 10.84, b = 352.8 and c = 1.794 ■ 10^. Here, the first term corresponds to the 
Madelung energy of the lattice, with = 0.895929 for a bcc-type lattice. The second 
term describes thermal contributions corresponding to an ideal gas of phonons while the 
higher-order terms account for anharmonic corrections. Assuming an unchanged ionic 
phase-space distribution function directly after the photoionization, where the spatial 
distribution due to the finite temperature is represented by Gaussian profiles centered 
on the sites of the bcc lattice, the internal energy of the initial nonequilibrium state is 
given by 



a2 




(4) 



where B = ^^a^, with B = 0.876719 for a bcc lattice. While the first two terms 
correspond to the ion kinetic energy and the Madelung energy of a bcc lattice of point 
charges, the last two terms give the leading-order corrections due to a broadened charge 
distribution at the lattice sites expanded in terms of the mean-squared ion displacement 
(r^). Assuming an initially non-interacting atomic gas and harmonic trapping potentials 
we have 



a2 Mula'^ ul 



where cul is the oscillation frequency of the trapping potentials. The final temperature 
is now simply obtained from energy conservation 

Tf 3/2 + u7i(Fi) 



Ti 3/2 + w,{V,) ' 



(6) 



with w = for F < 174 and w = for F > 174. Figure db) shows the ratio 
Tf/Ti as a function of Fi for (r^) /a^ = 0, i.e. for a point-like ionic density at the lattice 
sites. A comparison with our numerical data shows that equation © reproduces the 
numerical values for a finite-size plasma quite well. The remaining discrepancies should 
be attributable to finite-size effects |2I1I22]- At small values of Fi the temperature ratio 
starts out from unity and decreases with increasing Fj, since the excess potential energy 
of the final fluid state is naturally larger than that of the initial lattice configuration. 
As Fi exceeds a value of F; 70 the temperature ratio rises linearly, connected with a 
solidification of the plasma, where Ff = 174 stays constant (see inset). In this regime, 
regions of fluid and solid phases coexist in the plasma. At Fi ^ 84 the system has 
crossed the melting point and the temperature ratio approaches a value of 1/2. This 
factor of 1/2 originates from the additional thermal energy of 3/cB^f/2 arising from 
thermal lattice excitations. However, the situation changes if a finite extension of the 
initial charge distribution at the lattice sites is considered. Neglecting terms of order 
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Figure 2. (a) Final temperature as a function of the lattice filling factor for Fi = 200 
(solid line) and Fj = 400 (dashed line). The gray horizontal line separates cooling from 
heating, (b) Phase diagram for the final plasma state in the /-Fj plane. 



((r^) /a^) in equation (jH), one finds that for finite cul the temperature approaches a 
value of (1 + uji/ujl)/2 if (r^) < a\ 

Let us now turn to the question of the infiuence of possible lattice defects, i.e. how 
perfect must the initial-state lattice be in order to observe correlation cooling? Such 
defects are easily accounted for in the previous consideration. Introducing a filling factor 
/, defined as the probability that a given lattice site is occupied, simply reduces the 
average charge at each lattice site. Therefore, the potential energy of the initial state 
can be written as /^Wi(r*)/cBTi, where T* corresponds to the temperature scaled by 
the density of the perfectly filled lattice, which is related to the actual value of F; by 
r* = /"^/'^Fi. Substituting this in equation (jH) yields 

Tf ^ 3/2 + pw,if-'/^T,) 
Ti 3/2 + ^i;f(rf) 

as the generalization of equation (jHl) to non-perfect filling. In figure El^a), we show the 
resulting final temperature as a function of the filling factor for two different Fi. A 
decreasing filling enhances the initial potential energy (compared to that of a full lattice 
with the same average density, not the same lattice constant), leading to a rapid increase 
of the final temperature. As a consequence, cooling can be achieved for a nearly perfect 
lattice only. Furthermore, for a fixed / < 1, increasing Fi leads to larger values for 
the temperature ratio. This heating also infiuences the possibility of the creation of a 
crystallized plasma phase, since the melting point of Ff = 174 is shifted to larger values 
of Fj, as shown in the phase diagram figure Efb). For values of / < 0.9855, the system 
always ends up in the fiuid phase independent of Fj. 

The above picture suggests that, for a perfect initial state, cooling is always 
possible independent of the initial temperature. At low enough temperatures, however, 
the classical treatment will break down, and quantum mechanical effects will become 
important. As in the classical case the quantum statistical internal energy of a 
Coulomb solid|| can be written as a sum of the Madelung energy, harmonic contributions 

|j As mentioned previously, for typical experimental setups as in F^ is of the order of 10000, i.e. 
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arising from linear lattice excitations and anharmonic corrections. For the harmonic 
contributions to the excess energy we use an analytic approximation formula obtained 
recently |22j, which reproduces the numerical results for harmonic Coulomb crystals very 
accurately. The anharmonic corrections are taken from [21], where they were obtained 
by fitting the results of quantum Monte Carlo simulations. The resulting expression is 
too lengthy to be reproduced here, but incorporates both the known semiclassical high- 
temperature limit as well as the quantum mechanical groundstate energy of a Coulomb 
crystal, 

Wgs = -A^Tf + 1.329^^ - + O (4^] . (8) 

V^s \rs' ) 

Here, Q = (huo) / {k-QT) measures the importance of quantum effects, i.e. zero-point ion 
oscillations, and r^ = aMe'^/h'^ is the Wigner-Seitz radius in units of the ionic Bohr 
radius. As in the previous classical considerations, we apply a sudden approximation 
assuming an unchanged density matrix directly after the photoionization pulse. The 
initial energy of the plasma is thus found from the expectation value of the new 
Hamiltonian using a density matrix representing independent harmonic oscillators 
arranged on a bcc lattice. The result is 

2' 



a? 




(9) 



(10) 



which exactly corresponds to the classical expression equation (jH) except that the kinetic 
energy 3/cb71/2 is replaced by its quantum mechanical counterpart. Moreover, the mean- 
squared displacement is now given by 

rather than equation (0), where u = uq/ {uji^^yf\,). The final temperature is again 
obtained by equating the initial and final energies of the plasma. Figure 01 shows Ff 
as a function of Fi for r^ = 2 ■ 10^. In contrast to the classical case, here the final 
temperature cannot be decreased to arbitrarily low values, rather it approaches a finite 
constant value as Tj 0, since the initial charge distributions at the lattice sites have 
a finite extension at T; = due to zero-point oscillations. Using the low-temperature 
limit equation (jH)) as an estimate for the final energy and neglecting terms of order z/^ 
in the resulting energy balance, this limiting temperature is found to be 

^ = i^^ + l^-r], (11) 

where Q = 0.886 — 0.243/y¥^. Since Q < 0.886 independent of r^, the limiting 
temperature is always positive and takes on a minimum value of 

^^^^'^ =l-n'-^OAU (12) 



\ " / mm 

well within the crystallized regime. 
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Figure 3. Final scaled temperature as a function of the initial scaled temperature 
for rg — 2 ■ 10^ and for cjl/^o = 1 (thick solid line), ujl/^o = 2 (dashed line) and 
uji^/ujQ = 1/2 (dot-dashed line). For comparison, the result obtained with the classical 
expressions for initial- and final-state energies is shown by the thin solid line. The 
horizontal dotted lines show the estimate equation for the limiting temperature. 



Note, however, that the hmiting temperature obtained numerically with the full 
expression for W{ is slightly smaller since at finite temperature the system will not 
be in the groundstate, hence its potential energy is somewhat larger than that given 



by equation (jHJ. We may rewrite equation (fT^ as Ax/a = J2^T/[^yr^{l — il)] 



7.424r~^/^, where At = h/ \/2ti Mk-QT^ is the thermal deBroglie wavelength. This would 
suggest that it is possible to reach the degenerate regime of At/o ~ 1 if rg can be 
decreased sufficiently, i.e. to a value r^ ^ 3-10'^. Even for hydrogen this would correspond 
to an unrealistically large density of 4 ■ lO^^cm"'^. Note, however, that the value of r^ 
enters the ratio Xt/cl with a power of 1/4 only, so that values of Ax/a = 1/10 can 
already be reached if rg ^ 3 • 10'', corresponding to a density of 4 ■ lO^^cm"^, which 
has already been surpassed with cold Bose-condensed hydrogen atoms |2SI- In this 
regime, the present considerations will not hold, since effects of particle statistics will 
be important which have not been included in the above expressions for the energies. 
However, equation (fT^ might serve as an estimate for the parameter values necessary 
to reach the degenerate regime. 

In summary, we have discussed the possibility of correlation cooling and of the 
creation of a strongly coupled ultracold neutral plasma by photoionization of a highly 
ordered state of cold atoms, achievable with optical lattices. Classical considerations 
have shown that the initial temperature can be reduced by a factor of two if the atomic 
sample is noninteracting initially, offering the possibility to reach ion temperatures which 
have not been achieved so far. The quantum regime has been investigated using existing 
expressions for the equilibrium internal energy of the system. Here it was found that the 
final temperature approaches a constant value which arises from zero-point oscillations 



at an optimum ratio of cjl and of 




(13) 
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in the lattice potentials and depends on the density as well as the geometry of the optical 
lattice. An almost perfect lattice is necessary to observe the discussed cooling effects, 
since even a small amount of lattice defects largely increases the final temperature. 
However, it has been demonstrated that it is possible to create regular structures with 
a well-defined number of atoms localized in the potential wells of an optical lattice |2(ij . 

Financial support from the DFG through grant R01157/4 is gratefully 
acknowledged. 
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